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Abstract. We describe the application of the 'low-v' method to the extraction of the neutrino flux 
at MiniBooNE energies. As an example, we extract the relative energy dependence of the flux from 
published MiniBooNE quasielastic scattering cross sections with v < 0.2 GeV and V < 0.1 GeV 
(here V is the energy transfer to the target). We find that the flux extracted from the 'low-v' cross 
sections is consistent with the nominal flux used by MiniBooNE. We fit the MiniBooNE cross 
sections over the entire kinematic range to various parametrizations of the axial form factor. We 
find that if the overall normalization of the fit is allowed to float within the normalization errors, 
the extracted values of the axial vector mass are independent of the flux. Within the Fermi gas 
model, the distribution of the MiniBooNE data is described by a standard dipole form factor with 
Ma = 1.41 ±0.04 GeV. If nuclear transverse enhancement in the vector form factors is accounted 
for, the data are best fit with a modified dipole form factor with Ma = 1 . 10 =b 0.03 GeV. 
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PACS: 13.15.+g,13.60.-r,13.60.Hb,25.30.Fj,25.30.Pt 

(presented by A. Bodek at CIPANP 2012, St. Petersburg, FL, June 2012, and at NuFact 2012, 
Williamsburg, VA, July 2012). 

In a previous communication [1] we present the application of the 'low-v' method to 
the extraction of neutrino (v^) flux for energies {Ey) as low as 0.7 GeV. In this paper we 
extend the technique to Ey as low as 0.4 GeV and extract the relative energy dependence 
of the flux for the MiniBooNE experiment as an example. 

The charged current (v^) differential cross section can be written in terms of the 
square of the four momentum {Q^) and energy transfer (v) to the target nucleus. At low- 
V, if we integrate the cross section from v^/„ up to v= Vcut (where Vcut is small), we 
can write the 'low-v' cross section[l] in terms of an energy independent term which is 
proportional to the structure function W2, and small energy dependent corrections which 
are proportional to v/fi, or rrt^/E'^ where m — jiis the mass of the muon. 
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where 04 and 05 are negligible[l], and M is the nucleon mass. The uncertainties in the 
modehng of the small energy dependent correction terms are small, thus we can extract 
the relative flux from the number of low- v events at each Ey bin. 
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FIGURE 1. The ratio of the neutrino 'low-v' QE cross section (as a function of E^) to the 'low-v' QE 
cross section at = 1.1 GeV for v < 0.1 GeV (left) and v < 0.2 GeV (right). 
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FIGURE 2. The model uncertainty in the ratio of the neutrino 'low-v' QE cross section (as a function 
oiEy) to the 'low-v' QE cross section at = 1.1 GeV for v < 0.1 GeV (left) and v < 0.2 GeV (right). 

If we use the MINOS criteria that the fraction of events in the 'low-v' flux sample 
is lower than 60% of the total number of events in each E^y^ bin, we find that for 
neutrinos we can use events with v < 0.1 GeV to extract the relative flux for £'v > 0.4 
GeV, and events with v < 0.2 GeV for E^ > 0.7 GeV. For these v cuts, the cross section 
is dominated by quasielastic (QE) scattering. The flux extracted with the 'low-v' method 
is only a relative flux as a function of energy. It must be normalized at some energy. In 
this paper, we present the flux relative to the flux at £'v= 1.1 GeV. In our calculation of 
QE cross sections we use BBBA2007 electromagnetic form factors[3]. 

Figures 1 (v^) and 3 (v^) show the ratio of the 'low-v' QE cross section (as a function 
of Ey^y^) to the 'low-v' cross section at £'v,v^ = 1.1 GeV for v < 0.1 GeV (left) and 
V < 0.2 GeV (right) for various models. The data points are from the GENIE MC 
generator for a carbon target assuming a Fermi gas model and a dipole form for the 
axial form factor with Ma = 0.99 GeV[4]. The ratio is independent of the value Ma 
as illustrated by the fact that prediction of this ratio for a dipole axial vector mass Ma 
= 1.014 GeV (solid red line) and Ma = 1.3 GeV (solid blue line) are the same. Also 
shown are the changes in the prediction when we include nuclear enhancement in the 
transverse vector form factors [2] (TE) (shown as the solid black line). For £'v > 0.4 
GeV and v <0.1 and for Ey > 0.7 GeV and v <0.2 the ratio is approximately constant. 
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The uncertainties in the modehng of the energy dependence of this ratio are small as 
shown in figures 2 and 4. As an example, we extract the relative energy dependence of 
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fQ(1.1) for Antineutrino on Carbon v<0. 20 
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FIGURE 3. 

section at Ey 



The ratio of the 'low-v' QE cross section (as a function of Ey^) to the 'low-v' QE cross 
= 1.1 GeV for v < 0.1 GeV (left) and v < 0.2 GeV (right). 
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FIGURE 4. The model uncertainty in the ratio of the 'low-v' QE cross section (as a function of 
Ey^) to the 'low-v' QE cross section at E'y^ = 1.1 GeV for v < 0. GeV (left) and v < 0.2 GeV (right). 



the flux of the Fermilab booster neutrino beam from MiniBooNE data. The MiniBooNE 
experiment pubhshed[6] flux weighted double differential cross sections for QE neutrino 
scattering in bins of final state muon kinetic energy (E"^ — T^+m^) and muon angle 

(cos 0^). We extract the central value of v^^ = Ey^ — for each (r^, cos 0^) bin using 



2{M'^)E^-m'^f+ml-Ml) 



2^[{M'^)-E^ + jEl-mlco^e^] 



(1) 



where Mn and Mp are the neutron and proton mass, and M'^=Mn — Eb (Eb 
The left side of Fig. 5 shows the MiniBooNE bins of 0.1 GeV in 



= 34 MeV). 
and 0.1 in 



cos0^. The solid lines are lines of constant Ey 



QE 



The blue and red dotted lines are 
v^^ < 0.2 GeV, and v^^ < 0.1 GeV, respectively. We extract the 'low-v' flux from 
the MiniBooNE data as follows. Using the published MiniBooNE flux, we first fit the 
flux-weighted doubly differential cross section to three models. The parameters which 
are allowed to float within the models are the overall normalization and the axial vector 
mass Ma. The first model is a Fermi gas model with BBBA2007 electromagnetic form 
factors and a dipole form for the axial form factor. The second model includes Transverse 
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MiniBooNE 2D bins nu=o.2 




Muon Kinetic Energy T (GeV) Neutrino Energy (GeV) 

FIGURE 5. Left: The MiniBooNE QE cross section bins of 0. 1 in cos 0^ , and 0. 1 GeV in . The soHd 

Hnes are Hnes of constant Ey^. The blue and red dotted Hnes are v^^ < 0.2 GeV, and v^^ < 0.1 GeV, 
respectively. Right: The relative neutrino flux extracted from v^^ < 0.2 GeV cross sections (blue squares) 
and v^^ < 0.1 GeV (red squares) shown with statistical errors only. The black (flux A) and purple (flux 
B) lines are possible deviations (which are consistent with the 'low-v' flux) from the central values of the 
published flux. The green line is the quoted systematic uncertainty in the nominal MiniBoonNE flux. 
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FIGURE 6. Fa{Q^) measurements on free nucleons (a) Fa{Q^) re-extracted from neutrino-deuterium 
data divided by G^{Q^) (with Ma = 1.015 GeV). (b) Fa{Q^) from pion electroproduction divided by 
GoiQ^) corrected for for hadronic effects [7]. Solid black line - duality based modified dipole fit with 
Ma = 1.015G^y[3]. Short-dashed line - Fa{Q^)a2=v2- Dashed-dot line - constituent quark model[8]. Solid 
red line - duality based modified dipole with Ma = 1.10 GeV, which is our best fit to the MiniBooNE data 
on Carbon (accounting for Transverse Enhancement). 



Enhancement for the vector form factors[2] (BBC-TE) and a dipole form for the axial 
form factor. The third model includes TE for the vector form factors [2] (BBC-TE) and 
assumes a modified dipole form for the axial form factor as given in ref. [3]. The 
modification to the dipole form factor are from a fit[3] to all neutrino scattering data 
and pion electroproduction on free (H and D targets) nucleons as shown in fig. 6. The 
fit has the duality constraint that the vector and axial parts of structure function W2 for 
quasielastic scattering are equal at large Q^. 
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The ratio of the flux- weighted MiniBooNE measured cross sections at low-v to the 
calculated (with the nominal published MiniBooNE flux) flux-weighted cross sections 
for any of the three models is proportional to the ratio of the 'low-v' flux to the nominal 
MiniBooNE flux. As expected, the relative flux extracted as a function of neutrino energy 
is insensitive to the choice of model. 

The left side of Fig. 5 shows the ratio of the flux extracted from v < 0.1 GeV events 
(red circles) and v < 0.2 GeV events (blue squares) to the nominal flux. Only statistical 
errors are shown. The green line is the systematic error in nominal flux (as published by 
MiniBooNE). The extracted 'low-v' flux is consistent with the nominal flux within the 
quoted systematic errors. The black curve (flux A) and purple curve (flux B) are possible 
deviations (which are consistent with the 'low-v' flux) from the nominal flux. 

Next we fit for the best value of Ma for each of the three models. We find that if we 
let the overall normalization float within the systematic error the extracted values of Ma 
using the nominal flux, flux A, and flux B are within 0.015 GeV of each other as shown in 
Table 1. We find that with Transverse Enhancement, and a modified dipole form factor, 
the fit to the dependence of the MiniBooNE do /dQ^ on carbon favors an axial mass 
Ma = 1 . 10 ± 0.02 GeV. The ratio of this modified dipole fit with Ma = 1 . 10 GeV to the 
simple dipole parametrization with Ma = 1.015 GeV is shown as the solid red line in fig. 
6. The fit to the MiniBooNE data is more consistent with the values of Fa{Q^) extracted 
from pion electroproduction on free nucleons (shown in Fig. 6(b)), than with the values 
Fa{Q^) extracted from neutrino data on deuterium (Fig. 6(a)). 



TABLE 1. Fits to MiniBooNE neutrino quasielastic scattering data on carbon 



Form Factors 
vector/axial 


data set 
(2/D)/(lD) 


Ma 
(GeV) 


normalization 


model 


flux 


BBBA07 

FA=Dipole 


double diff (2D) 

dG/dQ^{m) 


1.35 ±0.02 
1.41 ±0.04 


0.99 ±0.01 
0.99 ±0.02 


39.8/135 = 0.30 
11.7/15 = 0.78 


nominal 
nominal 


BBC(TE) 
FA=dipole 


double diff (2D) 
dG/dQ^{m) 


1.22±0.02 
1.17±0.03 


1.01 ±0.01 
1.05 ±0.02 


43.6/135 = 0.32 
19.7/15 = 1.31 


nominal 
nominal 


BBC(TE) 
FA=mod. dipole 


double diff (2D) 
dG/dQ^{VD) 


1.17±0.02 
1.11±0.03 


1.01 ±0.01 
1.04 ±0.02 


35.2/135 = 0.28 
19.0/15 = 1.27 


nominal 
nominal 


BBC(TE) 
FA=mod. dipole 


double diff (2D) 

da/dQ^iVD) 


1.17±0.02 
1.10±0.03 


1.01 ±0.01 
1.04 ±0.02 


35.4/135 = 0.26 
17.8/15 = 1.18 


Flux A 
Flux A 


BBC(TE) 
FA=mod. dipole 


double diff (2D) 
do/dQ^{lD) 


1.17±0.02 
1.09 ±0.03 


1.01 ±0.01 
1.04 ±0.02 


38.3/135 = 0.28 
17.7/15 = 1.18 


FluxB 
FluxB 
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